Inhomogeneous states in two-dimensional frustrated phase separation. 
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We derive the phase diagram of a paradigmatic model of Coulomb frustrated phase separation 
in two-dimensional systems with negative short-range electronic compressibility. We consider the 
system subject either to the truly three-dimensional long-range Coulomb interaction (LRC) and to a 
two-dimensional LRC with logarithmic-like behavior. In both cases we find that the transition from 
the homogeneous phase to the inhomogeneous phase is generically first-order except for a critical 
point. Close to the critical point, inhomogeneities arrange in a triangular lattice with a subse- 
quent first-order topological transition to stripe-like objects by lowering the Coulomb frustration. 
A proliferation of inhomogeneities which have inside smaller inhomogeneities is expected near all 
the transition lines in systems embedded in the three-dimensional LRC alone. 
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I. INTRODUCTION 

Domain pattern formation is a beautiful example of 
cooperative behavior in complex systems with compet- 
ing interactions on different length scales^. In electronic 
systems, this idea has gained momentum due to theoret- 
ical and experimental studies in materials like cuprates 
and manganites. Indeed it has become clear that strong 
electron correlations generally produce a tendency to- 
wards phase separation in electron-rich and electron- 
poor regions^ii^iiii^iii^i^. Although the long-range part 
of the Coulomb interaction (LRC) spoils phase separa- 
tion (PS) as a thermodynamic phenomenon, the frus- 
trated tendency towards charge segregation is still im- 
portant and gives rise to inhomogeneous states where 
the charge is segregated over some characteristic distance 
but the average density is a fixed constant in order to 
avoid a diverging Coulomb cost in the thermodynamic 
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From the experimental point of view, a prominent 
emergent tool in recent years has been scanning tun- 
neling microscopy (STM). In the high temperature su- 
perconductors context, STM experiments in cuprates^^ 
have revealed a phase segregation on the scale of ^ 20 
lattice constants reconcilable with a Coulomb frustrated 
phase separation (FPS) between an underdoped pseudo- 
gap phase and a superconducting phase at higher dop- 
ings. Coexistence of insulating and metallic regions on 
the scale of tens to thousands of nanometers have been 
reported also in thin films of colossal magnetoresistance 
manganitesii. Noticeably percolation of the metallic re- 
gions is closely correlated to abrupt changes in trans- 
port suggesting that FPS is at the heart of the colossal 
magnetoresistance behavior^IiiS. More microscopic con- 
figurations with stripe patterns have been observed in 
cuprate&iS, nickelates^" and manganites^l. Even more, 
the discovery of pronounced anisotropies in transport 
measurements in ruthenates^S and in GaAs heterostruc- 



tures at weak magnetic fields^'^ are in line with the pro- 
posal of exotic electronic liquid phases^ analogue to 
the intermediate order phases of liquid crystalsS^. Ev- 
idence for mesoscopic electronic phase separation have 
been recently presented in organics2Si21, graphene^** and 
pnictides^^. 

The generality of this phenomena calls for simple mod- 
els which neglect the microscopic details of the specific 
system rather capturing its general properties. The ten- 
dency towards PS is then caused by the appearance of 
anomalies in the electronic contribution to the free energy 
of the system /e, as a function of the density n. Two kind 
of anomalies have been identified as the most relevant 
ones for phase separation in electronic systemsiiii^i^. 
The first anomaly consists of a range of densities where 
the compressibility is negative while the second one cor- 
responds to a single point where the inverse of the elec- 
tronic compressibility has a Dirac-delta-like negative di- 
vergence due to a crossing of the free energies of the two 
competing phasesi^ii^i^. The two anomalies can be la- 
beled by the exponent 7 characterizing the behavior of 
the electronic free energy around a reference density ric, 
i.e. fe = OL\n — nc\'^ with a < and 7 = 2 (negative 
compressibility region) or 7 = 1 (cusp behavior). Both 
anomalies appear often in model computations of elec- 
tronic systems, for example Refs. IsHol flolls^ correspond to 
7 = 2 while Refs. i33ii34 to 7 = 1. 

When LRC effects can be considered as a weak pertur- 
bation upon the PS mechanism, one can achieve a uni- 
versal picture of the FPS phenomenon. Antithetically, 
upon strong frustrating effects, the two short-range com- 
pressibility anomalies give rise to two different universal- 
ity classes. Within the 7 = 1 universality class, differ- 
ent studies in 3D systemsiii^ and in 2D systemsi^ii^i^ 
have enlightened the key role of the system dimension- 
ality in the FPS mechanism. Recently, the full phase 
diagram for the 7 = 2 class has been characterized in 
Ref. |3^ both in isotropic and strongly anisotropic 3D 
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systems. The aim of this work is to investigate the prop- 
erties of the FPS phase diagram for 7 = 2 at strong 
frustration in 2D systems. This is an important gener- 
ahzation as many interesting strongly correlated systems 
are layered as cuprates, nickelates, iron oxypnictides su- 
perconductors, or are truly two dimensional as the two- 
dimensional electron gas in heterostructures, graphene, 
etc. Relevant for cuprates, manganites and nickelates is 
how charge density wave instabilities of a uniform phased 
evolve into the very anharmonic structures experimen- 
tally observedi2i2£i2i. 

We will consider two versions of the Coulomb interac- 
tion labeled by an effective dimensionality d. The d = 3 
case corresponds to the usual three dimensional Coulomb 
interaction and applies to the physical systems above 
mentioned (layered systems, heterostructures, graphene, 
etc.). The d = 2 case, instead, corresponds to a fictitious 
logarithmic Coulomb interaction.-^ The effective dimen- 
sionality d regards only the interaction and should not 
be confused with the dimensionality of the system which 
is 2D. Alternatively the 2D system with logarithmic in- 
teraction can be considered as an anisotropic 3D system 
subject to the conventional Coulomb interaction where 
modulations in one "hard" direction are forbidden. This 
anisotropy can originate in the underling crystal struc- 
ture and the case of only one hard direction (as opposed 
to two hard directions considered in Ref . f36l ) corresponds 
to the case of weakly coupled chains. This is relevant for 
electronic phase separation in organics where elongated 
domains have been reported i^^'^^ 



II. PHASE DIAGRAMS IN 2D SYSTEMS 

Our starting point is a paradigmatic model of Coulomb 
frustrated phase separation in systems with a short-range 
negative compressibility density region that is defined by 
the following free energy: 

F = J dx c\Vn{x)\^ + a[n{x.) — Ticf + l3 [n(x.) — Tic]'^ 
I [n(x) — n] [7i(x') — n] 



-I- — I dn I dx 



G-i(lx-x'l) 



(1) 



Here the gradient term models the surface energy of 
smooth interfaces and is parametrized by the stiffness 
constant c whereas a < is proportional to the inverse 
short-range compressibility. The inclusion of a fourth or- 
der term with /3 > provides a symmetric double-well 
form of the short-range part of the free energy with min- 
ima at n = ric ± [|q;|/(2/3)]^/^. In general, a will depend 
on external parameters like the pressure. It can be also 
taken as temperature-dependent, as in Landau theory, in 
which case the model becomes a mean-field description 
of a temperature-driven transition to an inhomogeneous 
state. In addition, Sq indicates a static dielectric constant 
due to external degrees of freedom and n is the average 
density. A rigid background ensures charge neutrality. 



Finally G(|x — x'|) is a positive-definite LRC kernel that 
in the d — 2 case has a logarithmic-like behavior whereas 
for d = 3, G~i(|x - x'l) = |x - x'|. 

Since the model Eq. ^ has several parameters, it is 
convenient, in order to study the phase diagram, to mea- 
sure all lengths in unit of the bare correlation length 
^ = ^/2c/a, densities in units of [a/ {2(3)]^/^ and energy 
densities in terms of the barrier height /o = a^/(4/3). 
Then one reaches a dimensionless form consisting, apart 
from an irrelevant constant, of a model augmented 
with a long-range Coulomb interaction. The correspond- 
ing Hamiltonian is defined by: 

n = J rfx[02(x)-l]%|V0(x)|2 + ^ 

dx I .x'm=am-l^ (2) 



G-i(lx-x'l) 

where, is dimensionless renormalized frustration pa- 
rameter of the long-range interaction given by 
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and the classical scalar field (f) represents the dimension- 
less local charge density with average density 4>. The 
effect of frustration can also be measured_by the param- 
eter introduced in Refe. 35,1223. It can be 
obtained from the following dimensional analysis. Tak- 
ing (/) ~ 1, the Coulomb energy density per domain 
can be estimated making the integrals in a volume of 
order as Q^l'^~'^ with Id indicating the typical size 
of the inhomogeneities measured in unit of ^. Within 
a sharp interfaces approach, the surface energy density 
goes like 1 1 Id- Both quantities are optimized whenever 
Id ~ l/Q^/^^-'^\ Then g2/(5-d) ^^i^gg ^Ylg meaning of 
ratio between the additional energy cost induced by frus- 
tration and the typical PS energy gain. 

For both the truly long-range Coulomb interaction and 
the logarithmic-like interaction, the charge susceptibility 
in momentum space can be derived by computing the 
static response to an external field. At k ^ one finds: 



X(k,rf) = 



i4-d 



-2 



(3) 



Notice that the second term in the brackets yields the 
familiar forms of the Fourier transform of the Coulomb 
interaction for d = 2, 3. 

The charge susceptibility has a peak at finite momen- 
tum fco = [(4 - d)22-'^7rg2] ^'^^-^^ determined by the 
competition between interface and charging effects. It 
diverges as Q approaches a Gaussian instability line : 
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geneous phase can be expressed in momentum space as: 
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FIG. 1: (Color online) (a) Contour plot of the charge density 
modulation [4>{'x.) — 0] /(t>G near the critical point for a trian- 
gular crystal of inhomogeneities. The yellow spots are in the 
vicinity of the triangular lattice points where ^(x) = <p + 6(j)G- 
(b) Same for a unidirectional SCDW. In both panels, lengths 
have been measured in units of 1/ko for convenience. 



represents the maximum frustration degree above which 
the system is in its homogeneous phase. 

The Gaussain instability Une Eq. ([4]) indicates an insta- 
bility towards a sinusoidal charge density wave (SCDW) 
of period 27r/fco and direction chosen by spontaneous 
symmetry breaking. Still, it cannot persist up to Q — > 
since inhomogeneities are predicted only for average den- 
sities inside the "miscibility" gap^S I*/*! < l/v^- This is 
not in tune with the usual Maxwell construction that on 
the contrary predicts an inhomogeneous state for |0| < 1. 
Analogously to 3D systems'^^ indeed, we find that the in- 
clusion of non-gaussian terms results in a first-order tran- 
sition preempting the Gaussian instability except for the 
critical point (CP) (0, Q) — (0, Qc) as we show below. 

By restricting to periodic textures, the free energy den- 
sity difference between a modulated phase and the homo- 
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— = 5] 0GX"l(G,d)0_G 



(5) 



G#0 



Gi,G2,G35^0 

+ ^ </'Gi0G2'/'G3'?!'G4^Gi + G2+G3+G4,O 

Gi,G2,G3.G4#0 

where the G's are the wavevectors of the reciprocal lat- 
tice and V indicates its unit cell volume. The presence 
of the cubic term in Eq. ^ opens the way to first-order 
transitions except for the CP. Away but close to the CP, 
the transition will be weakly first-order. Hence it can 
be treated, analogously to the liquid-solid transitio n^^'^^ , 
by fixing the wavevectors to have an equal magnitude 
|G| = fco. In this case, to gain an energetic advan- 
tage from the cubic term, one needs triads of wavevec- 
tors forming an equilateral triangle so that the Kronecker 
delta is satisfied. In 2D systems this is only verified by 
the hexagonal reciprocal lattice that is defined by the set 
of six wavevectors G/ fco = (±cos7r/3,±sin7r/3), (±1,0). 
The corresponding free energy density reads: 
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(6) 



As a result, one finds that the first structure to become 
stable corresponds to droplet-like structures forming a 
triangular crystal of inhomogeneities. Fig. [T] (a) shows 
the corresponding charge modulation close to the critical 
point. 

Upon minimizing Eq. ([5]) with respect to the wave am- 
plitude (j}Q and requiring 5F /V = 0, one finds that the 
first-order transition line is given by: 
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On entering in the inhomogeneous phase density re- 
gion, the triangular crystal phase is expected to compete 
with a unidirectional SCDW in which the local charge is 
modulated along stripes [see Fig.[ljb)] whose free energy 
density reads: 



5F_ 
V 



20G 



64 



(7) 



We find that at fixed frustration Q, stripes become 
stable close to cf> = [see Fig. [2] . This leads to a first- 
order morphological transition that partially restores the 
translation symmetry. In the frustration-density plane, 
the morphological transition line can be determined by 
equating the free energies densities Eqs. ([21), ([7]) of the 
two phases minimized with respect to the wave amplitude 
0G- We find : 
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FIG. 2: (Color online) Comparison for Q — 0.9 Qc of the free 
energy densities of the triangular crystal phase (full lines) 
and the SCDW (dashed hues) for a d = 2 LRC (thick lines) 
and the truly d = 3 LRC (thin lines). We have substracted 
the free energy density of the homogeneous phase (horizontal 
line). Arrows indicate the cusp singularity induced in both 
cases at the topological transition. 

Although the transition hnes are only asymptotically 
exact close to CP—, we find the topology of the phase 
diagram found in the strong frustration regime to per- 
sist at weak frustration. This is shown in Fig. [3] for 
a d = 2 Coulomb interaction where we characterized 
the full phase diagram. In the weak frustration regime 
{Q << 1), we considered a Uniform Density Approxi- 
mation (UDA) where we assumed domains of uniform 
density of one or the other phases separated by sharp 
interfacesiiii^i^^'^i^S. The UDA is an accurate descrip- 
tion in this regime since there is a sufficient separation 
among the typical size of the domains Id and the screen- 
ing length Is that controls the charge relaxation inside the 
domainsii. For d = 2,3, indeed, the latter goes like--: 



In addition, is much larger then the typical Gaussian 
instability wavelength 1/fco. Thus one can rely on a sharp 
interface approach where surface energy effects are taken 
into account by neglecting long-range effects and com- 
puting the excess energy of an isolated interface. Then 
one recovers the surface tension of a kink^^ a = 8/3. 

In the weak frustration regime, we find the same topo- 
logical transitions (thick lines in Fig. [3]) as in strong frus- 
tration but now the inhomogeneities form sharply defined 
circular droplets and stripes. We have also studied the 
crossover from weak to strong frustration minimizing a 
discretized version of the model Eq. ^ in the Wigner- 
Seitz Approximation. The corresponding transition lines 
are shown with squares and diamonds in Fig. [3l For 
droplet-like structure we assumed circular symmetry in 
order to reduce the effective dimensinality of the mini- 
mization procedure to one. The upper panel of Fig. [4] 
shows the corresponding charge density profile in the 
radial direction at = 0.5 and different frustrations. 
By decreasing the frustration degree, the charge density 
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FIG. 3: (Color online) Phase diagram in two-dimensional sys- 
tems embedded in the d = 2 long-range interaction. The 
small dots indicates the Gaussian instability line Qg. The 
thin (thick) lines represent first-order transitions in the strong 
(weak) frustration approximation. In the two limits they over- 
lap with the corresponding numerically determined transition 
lines from the homogeneous phase to droplet inhomogeneities 
(□) and from droplets to stripes (<0). Finally the black circle 
indicates the CP. 



modulation in the radial direction gets more unharmonic 
similarly to what we find for striped states at cj) — [see 
lower panel of Fig. 2]. 

The features of the FPS phase diagram are dramati- 
cally different if the stiffness constant is made strongly 
anisotropic. More in detail, we assume that c\\/c± > 1 
where C|| {c±) indicates the stiffness component in the 
"hard" ("soft") direction of the 2D plane and focus on 
the limit where any charge modulation in the "hard" di- 
rection is completely forbidden, i.e. C||/cj^ oo. In a 2D 
system subject to the logarithmic-like d — 2 interaction, 
in view of the analogy discussed in Sec. IH this would cor- 
respond to consider a 3D system with two "hard" direc- 
tions. The corresponding phase diagram has been shown 
in Fig. 3 of Ref. [Sy. In that case it has been shown that 
this strong anisotropy allows for second- and first-order 
transitions transition lines joined by a tricritical point. 

We now show that this feature appears also in 2D sys- 
tems subject to the ordinary d = 3 Coulomb interaction. 
Since the strong anisotropy allows only for unidirectional 
charge modulations, the cubic term of Eq. ^ identically 
vanishes by restricting to wavevectors of feg magnitude. 
Thus, to include non Gaussian terms, one has to retain 
at least two collinear wavevectors so that the Fourier de- 
composition of the order parameter reads: 

0(x) = (/) + 201 cos (Gi • x) + 202 cos (G2 • x) 

where IG2I — 2|Gi| = 2fco. By minimizing with respect 
to the periodicity of the charge density modulation and 
assuming 02 << 0i, that can be checked a posteriori, 
one obtains the following Landau free energy expansion 
in terms of 0i: 

—^ = r^l+ Ui4>\ + "60? (8) 
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FIG. 4: (Color online). Top panel: Behavior of the lo- 
cal charge density modulation for droplets at — 0.5 and 
different values of the Coulomb frustration Q/Qc- As the 
Coulomb frustration decreases, unharmonicity is build in and 
the charge density profile tends to a square wave in the radial 
direction r (thick full line). Bottom panel: Same for <^ = 
where only stripe-like objects appear. 



where the quadratic term coefficient 



r = 4 



Q_\ 

QcJ 



4/3 



Qc 
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vanishes along the Gaussian instabihty line whereas U4 = 

6 — 540 {Qc/QY^^ and uq is a positive constant, ft 
is simple to notice that this is the canonical form of 
the Landau free energy expansion around a tricritical 
point that is determined by the vanishing of U4 along the 
Gaussian instability^^ and is therefore given by {</), Q) — 
(1/(2V3), (3/4)3/4). rpj^g density window of the 

second order line is important because it allows for a 
second order quantum critical point (QCP) as discussed 
m Ref's'. Furthermore as parameters are changed (den- 
sity, frustration) an evolution similar to the one displayed 
in Fig. |4] will occur. This establishes a connection be- 
tween the Gaussian instabilities predicted in models of 
the cupratesfi and the probably anharmonic stripe pat- 
terns observedi^. 

Apart from small shifts of the transition lines, the ef- 
fective dimensionality of the Coulomb interaction does 
not change the features of the 2D phase diagram both 
in absence and in presence of strong anisotropics. This 
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FIG. 5: (Color online). The phase diagram for anisotropic 2D 
systems subject to the truly d — 3 Coulomb interaction. The 
thin line corresponds to the Gaussian instability that deter- 
mines the second-order transition line above the tricritical 
point (black circle). Below the tricritical point, the transition 
is first-order and is determined in the weak coupling regime 
assuming the UDA (dashed lines). The long-short dashed 
lines are qualitative lines of the crossover from the weak to 
the strong coupling regime. 



is reported in Fig. [5] where we show the phase diagram 
of an anisotropic 2D system subject to the truly d = 3 
Coulomb interaction. An important difference, however, 
will appear close to the transition lines as discussed be- 
low. 



III. DISCUSSION AND CONCLUSIONS 

The topology of the FPS phase diagrams in 2D systems 
has a strong similarity to the 3D case of Ref. HH except 
that a crystal of drops appears in addition in 3D. Thus 
one can state that the system dimensionality D and the 
effective dimensionality of the long-range Coulomb inter- 
action d do not qualitatively affect the FPS phenomenol- 
ogy in the 7 = 2 universality class. This is not quite the 
end of the story since the transition from the homoge- 
neous to the inhomogeneous phase and the topological 
transitions are first-order-like. A cusp singularity in the 
minimal free energy of the system will be naturally pro- 
duced as indicated by arrows in Fig. [21 This drives the 
system towards a new stage of FPS tendency which now, 
however, is governed by the 7 = 1 universality class and a 
new renormalized frustration constant Q' which depends 
on the distance from the critical point and diverges at 
the critical point. At this second stage of domain pat- 
tern formation, a dramatic difference between 3D and 
2D d — 3 systems will arise. Indeed in 3D systems, for 
7 = 1, there is a critical value of Q' above which the 
second stage of pattern formation is suppresse d^^i^^'^'^i'^^ 
while in a 2D, d — 3 system the second stage always 
occurs no matter how big is Q' For the case of 2D 
d ~ 2 systems, a critical frustration for Q' exists since, 
as mentioned in Sec. HI it can be viewed as an anisotropic 
3D, c? = 3 system. 
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The stabilization of inhomogeneities which have inside 
smaller inhomogeneities represents a relevant new effect 
to be considered in the 2D, 7 = 2 phase diagram. These 
states will consist of stripes of droplets alternating with 
stripes of stripes (homogeneous stripes) and will appear 
in a narrow but still finite (j) window near the first-order 
morphological (homogeneous-inhomogeneous) transition. 

In conclusion, we have derived the Coulomb frustrated 
phase separation phase diagram in 2D systems with a 
short-range negative compressibility subject either to a 
fictitious two-dimensional logarithmic-like Coulomb in- 
teraction and the truly three-dimensional Coulomb inter- 
action. Similarly to three-dimensional systems, we find 
that the transition from the homogeneous phase to the 
inhomogeneous phase is always first-order except for a 
CP. Close to the CP, inhomogeneities are predicted to 



form a triangular lattice with a subsequent transition to 
striped states. The transition lines continuously evolve 
into the weak frustration limit. In 2D systems subject 
to the three-dimensional Coulomb interaction both first 
order and second order transition lines are found sepa- 
rated by a tricritical point. A proliferation of inhomo- 
geneities is naturally expected near the first order transi- 
tion lines. The second order transition lines make a QCP 
when crossed by the change of one parameter, like doping 
considered in Ref. [^. Then one can expect very harmonic 
charge modulations close to the QCP which evolve into 
more anharmonic modulations as the distance from the 
QCP increases. 
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